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ABSTRACT 

We present long-slit spectroscopic data which reveals extended Ly a emission around the z = 6.417 radio- 
quiet quasar CFHQS J2329-0301. The Ly a emission is extended over 15 kpc and has a luminosity of > 8 x 
10 36 W, comparable to the most luminous Ly a halos known. The emission has complex kinematics, in part due 
to foreground absorption which only partly covers the extended nebula. The velocity ranges from -500 km s -1 
to +500 km s" 1 , with a peak remarkably close to the systemic velocity identified by broad Mg II emission of 
the quasar. There is no evidence for infall or outflow of the halo gas. We speculate that the Lya emission 
mechanism is recombination after quasar photo-ionization of gas sitting within a high-mass dark matter halo. 
The immense Ly a luminosity indicates a higher covering factor of cold gas compared to typical radio-quiet 
quasars at lower redshift. 

Subject headings: cosmology: observations — quasars: general — quasars: emission lines 
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1. INTRODUCTION 

There are still many mysteries surrounding the cosmologi- 
cal evolution of the most massive galaxies. Did these galaxies 
undergo a monolithic collapse with extremely high star for- 
mation efficiency at an early epoch? Did they all build up 
supermassive black holes rapidly at their centers? What pro- 
cesses are responsible for the subsequent halt in star forma- 
tion, so that today they appear red and dead? 

At high redshift, massive galaxies are extremely rare and 
are not found by small area surveys such as with the Hub- 
ble Space Telescope. One of the best ways to locate mas- 
sive galaxies at z ~ 6 is by finding luminous quasars, since 
these quasars contain supermassive black holes with masses 
> 10 8 M Q accreting at the Eddington limit (Jiang et al. 2007; 
Kurk et al. 2007; Kurk et al. 2009; Willott et al. 2010). The 
space density of black holes at this epoch is a factor of 10 4 
times lower than today (Willott et al. 2010). They are located 
in rare high matter density peaks where galaxy formation got 
underway rapidly (Volonteri & Rees 2006). 

Cosmological simulations show that massive galaxies lo- 
cated within hot gas halos can accrete a substantial amount 
of cold gas, so long as the cold gas is fed into the galaxy in 
narrow filaments (Dekel et al. 2009). These cold streams 
may be the dominant mode of star formation at high redshift. 
This is difficult to prove observationally by Ly a absorption, 
due to the small covering factor of the cold streams (Faucher- 
Giguere & Keres 201 1). 

Cold gas accretion can also have effects on the Ly a emis- 
sion around galaxies. As cold, neutral gas falls toward the 
central galaxy it can be heated and ionized by shocks and re- 
combines emitting Ly a photons (Haiman et al. 2000). Cold 
gas within galaxies could also be ionized by an intense star- 
burst or an active galactic nucleus (AGN) at the galaxy center 
(Haiman & Rees 2001). Ly a "blobs" or halos have been dis- 
covered around active and inactive galaxies (e.g. Heckman et 
al. 1991; Fynbo et al. 1999; Steidel et al. 2000; Matsuda 
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et al. 2004; Weidinger et al. 2004; Christensen et al. 2006; 
Hennawi et al. 2009). All of the above mechanisms have been 
proposed to explain the Ly a emission from these objects and 
it is likely that a range of mechanisms are at work. 

CFHQS J232908-030158 (hereafter CFHQS J2329-0301) 
is a quasar at z = 6.417 discovered in the Canada-France High- 
z Quasar Survey (CFHQS; Willott et al. 2007). It is radio- 
quiet (Wang et al. 2008) and powered by a black hole with 
mass w 2 x 10 8 M Q (Willott et al. 2010). There is some ev- 
idence for a protocluster of companion star-forming galax- 
ies based on deep multi-color imaging (Utsumi et al. 2010). 
However, most of these objects are visible in the i' band 
thumbnail images and therefore likely lie at redshift lower 
than z = 6.4. Goto et al. (2009; hereafter G09) presented deep 
broad-band imaging observations of this quasar and showed 
that it is spatially extended at z 1 band and possibly also at z r 
band. The high level of extended emission at z' band sug- 
gests that there is a large flux due to extended Ly a emission. 
Willott et al. (2007) found that the broad Ly a line in this 
quasar has twice the equivalent width of typical quasars. 

We present a deep long-slit spectroscopic observation of 
CFHQS J2329-0301 with the aim of determining how much 
of the extended flux of G09 is due to Ly a emission and its 
spatial and kinematic distributions. In Section 2 we analyze 
the imaging data of G09. In Section 3 we describe the spectro- 
scopic observations and their analysis. Section 4 discusses the 
physical nature of this galaxy. All optical and near-IR magni- 
tudes in this paper are on the AB system. Cosmological pa- 
rameters of H = 70 km s" 1 Mpc" 1 , n M = 0.28 and fl A = 0.72 
(Komatsu et al. 2009) are assumed throughout. One arcsec- 
ond on the sky corresponds to a physical size of 5.6 kpc at 
z = 6.417. 

2. EXTENDED EMISSION IN BROAD-BAND IMAGES 



G09 presented 
CFHQS J2329-0301 
and z r - The data 
with seeing of 0.5". 



Subaru Suprime-Cam images of 
in three broad-band filters: /', z' 
were obtained in excellent conditions 
The quasar is very faint in the i' band 
(/' = 25.5) due to IGM absorption shortward of the Lya 
emission line. Therefore no extended emission in this filter 
is observed. However, the data are useful to check that none 
of the extended emission found at longer wavelengths is due 
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FIG. 1 . — ie/f-' Subaru Suprime-Cam z' band image of the quasar CFHQS J2329-0301 from G09. The red lines marks the location of the Keck ESI 1" wide 
spectroscopic slit. Middle: The same image after subtraction of a model PSF scaled to 0.87 times the total flux. Right: The PSF-subtracted image after smoothing 
with a 1 pixel (0.2") gaussian kernel. All three images are displayed with the same greyscale stretch parameters. 



to a contaminating object such as a low redshift galaxy. At 
zf band the object is clearly more resolved than stars in the 
field. G09 carried out point-spread function (PSF) subtraction 
revealing an elongated (roughly E-W) object extended over 
4 arcsec with magnitude z' = 23.5 ±0.3 in a 2.6" radius 
aperture. Due to the lower sensitivity of Suprime-Cam at z r 
band, G09 found only a tentative detection of extended flux 
in this filter (which probes longward of the Lya emission 
line) with magnitude z r > 24. The z' band contains the Ly a 
line and the brightness of the extended emission in this filter 
suggests that at least some fraction of the extended emission 
is from Ly a. 

In order to compare these imaging observations with our 
spectroscopy to be presented in the following section, we have 
obtained the processed Suprime-Cam z' band image from the 
authors of G09. We performed an independent PSF subtrac- 
tion using field stars of similar magnitude to the quasar. Six 
stars were used to construct a model PSF which has FWHM 
2.4 pixels or 0.48". The PSF was scaled to match the total 
magnitude of the quasar in a circular aperture of radius 2". 
This PSF was then subtracted from the quasar image. By trial 
and error it was found that the model PSF should be scaled by 
a factor of 0.87 such that the integrated flux within the central 
few pixels was just above zero. Note that this method (maxi- 
mal PSF-subtraction) therefore gives a lower limit on the ex- 
tended flux. We determined an upper limit to the PSF scaling 
by considering when the integrated flux within the central few 
pixels was more negative than the noise across clean parts of 
the image. It was found that the PSF scaling had to be < 0.89 
at the 4 a level. 

The PSF-subtracted residual image (center panel of Figure 
[T]) shows a similar structure to that in G09. The ring-like na- 
ture of the residuals could be real or could be an artifact of 
the maximal PSF-subtraction. We measure a magnitude of 
z' = 23.45 ± 0.06 within a 2" radius aperture, consistent with 
the magnitude measured by G09. Note the uncertainty on the 
magnitude quoted above is only due to the measured pixel- 
to-pixel noise in the sky aperture. It does not include the un- 
certainty due to PSF scaling or the effect of correlated pix- 
els. The 4 a limit based on the PSF scaling discussed previ- 
ously corresponds to an extended component magnitude limit 
ofz'< 23.62. 

Inspection of the image suggests substantial correlation be- 
tween nearby pixels due to a spatially variable background. 



To ascertain the significance of the extended emission in the 
Z 1 band image, we placed 2" radius apertures at random across 
the image (after masking of real sources) and measured the 
distribution of fluxes in these apertures. A flux equivalent to 
magnitude fainter than zf = 25.2 is found in 68% of the ran- 
dom apertures and magnitude fainter than zf = 24.6 is found 
95% of the time. Assuming a Gaussian distribution, we ex- 
trapolate this to infer that the extended z' emission has a sig- 
nificance of 5 a. G09 stated that the extended z' emission has 
significance 16 a. However, their random aperture results are 
similar to ours; they found a la sky noise of z' = 24.9, which 
would make the detection significant at only 3.6 a. In con- 
clusion, the extended emission at z' band is significant at at 
least the 5 a level (recalling that we subtracted off a maximal 
PSF). With a magnitude of zf = 23.5 the extended emission 
comprises at least 13% of the total flux from this object in this 
filter passband. 

3. LONG-SLIT SPECTROSCOPIC OBSERVATIONS OF THE 
EXTENDED EMISSION 

3.1. Observations 

The quasar CFHQS J2329-0301 was observed with the 
Echelle Spectrograph and Imager (ESI) spectrograph (Shei- 
nis et al. 2002) at the Keck-II telescope during the nights of 
4+5 October 2007 and 24 September 2008. The atmospheric 
conditions were variable and only the data taken in reasonable 
conditions were included. The total integration time was 8.5 
hours, split into 30 minute exposures with the quasar position 
offset along the slit between each exposure. The 1" wide slit 
was positioned at an angle of 105° East of North, as indicated 
on Figure [T] Note that the spectroscopic observations were 
performed before publication of G09, so the good alignment 
of the slit with the elongation angle of the extended emission 
is coincidental. Ly a falls in order 7 of the echelle, where the 
spatial pixel scale is 0.163" pixel" 1 and the spectral scale is 
11.5 km s" 1 pixel" 1 . The spectral resolution with the 1" slit is 
75 km s" 1 or fl = 4000. 

Data reduction was performed using the ESIRedux codeQ 
After bias removal, a pinhole mask observation was used to 
determine the spectral traces of the orders on the detector. The 
data were then flat-fielded, sky-subtracted and the wavelength 
solution determined using arcs. A standard star observation 

5 http://www2.keck.hawaii.edu/inst/esi/ESIRedux/index.html 
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FIG. 2.— f/pper: Keck ESI extracted spectrum of CFHQS J2329-0301. The 
velocity axis is for the Ly a line relative to the systemic redshift of z = 6.417. 
Middle: The spatial Gaussian FWHM fitted to each column of the spectrum 
in units of kpc and arcsec. Regions with S/N < 10 in the extracted spectrum 
are marked with grey shading, indicating difficulty in fitting the Gaussian. 
The dashed line show the best fit linear solution to good wavelength regions 
as described in the text. The dotted line at velocity -190kms~' marks the 
location of an absorption feature in the broad Ly a line with associated N V 
absorption. At the velocity of this Ly a absorption feature the flux is more 
spatially extended than the continuum and broad Ly a emission. Lower: Spa- 
tial location along the slit of the peak of the fitted Gaussian of each column. 
The dashed line is the best fit linear solution to good wavelength regions. The 
spatial centroid of the flux at velocity -190 kms~' is offset by 0.2" (~ 1 kpc) 
from the continuum and broad Ly a. 

was used to determine the target traces. There is significant 
fringing at the red end of the spectra, so a master fringe frame 
was produced by median-combining all data of faint objects 
obtained over the two runs with the locations of known tar- 
gets masked out. Appropriately scaled fringe frames were 
subtracted from each spectrum of the quasar. The individ- 
ual frames were then shifted and combined with weighting 
dependent upon their signal-to-noise ratios (S/N). Each order 
was extracted from the 2D image and rectified using the target 
trace. The seeing, as measured from the quasar continuum is 
0.93", which is significantly poorer than in the imaging ob- 
servations of G09. Flux calibration is achieved with standard 
star spectra and the published broad-band magnitude of the 
quasar (Willott et al. 2007). 

3.2. Spectral PSF subtraction 

Due to the relatively short slit of ESI (20") and the need 
to dither the target along the slit, it was not possible to si- 
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FIG. 3. — Similar series of plots as in Figure f2]but focusing on the region 
around the hya line. Uncertainties on the fitted Gaussian FWHM in the 
middle panel are plotted. The absorption feature at velocity -190kms -1 is 
more clearly revealed to have a complex structure. The broad Ly a line has a 
FWHM that is about 5% greater than the linear fit to the continuum. 



multaneously observe a star on the slit to determine the PSF. 
Therefore the true PSF for these observations is uncertain. In- 
spection of the spectrum showed that there is extended flux 
at the wavelength of the Ly a line. Figures [2] and [3] show 
the extracted one-dimensional spectrum and how the spatial 
FWHM and peak location along the slit vary as a function of 
Ly a velocity. Zero velocity is set at z = 6.417, based on the 
broad Mg II emission line redshift from Willott et al. (2010). 
The Mg II line is a good estimator of the systemic velocities 
of quasars (Richards et al. 2002). 

There are several deep Lya absorption lines on the blue 
wing of the broad Ly a line. The broadest, deepest H I ab- 
sorption dip centred at -190km s" 1 (z a b s = 6.4123) is a multi- 
component H I system. Fitting the Ly a profile yields a mini- 
mum of three (more likely four) components, of which one 
has associated Nv absorption. The total neutral hydrogen 
column density is logA^(Hl)= 14.8 ± 0.2cm" 2 . The asso- 
ciated Nv doublet at z abs = 6.41241 ±0.00005, marked on 
Figure |2j is unresolved and highly saturated (both lines of 
the doublet have rest equivalent widths w r = 0.20 A). There 
is also possible associated Si II 1260.4 absorption at z a bs = 
6.4136 ± 0.0001 but it falls in a region of strong sky absorp- 
tion lines. 

A Gaussian flux distribution was fitted to each column of 
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the spectrum to determine the spatial FWHM and peak loca- 
tion (centroid). It was found that both these quantities were 
slowly varying functions of wavelength. Linear fits to both 
quantities were determined using a % 2 minimisation using 
only wavelength regions with high quality data and avoid- 
ing the obviously more spatially extended parts of the spec- 
trum. The fitted regions were -700 to -1700 km s -1 and +700 
to +8440 km s" 1 . 

Figures |2] and |3]reveal two interesting things. Firstly, at the 
location of the H I absorption at -190 km s" 1 , the Ly a flux is 
significantly more extended than elsewhere with a FWHM of 
1.2 arcsec. Secondly, its centroid is significantly offset along 
the slit by 0. 1" to 0.2" (« 1 pixel) towards the west. Note that 
this offset is relatively small compared to the total extent of 
extended emission found in the image of G09. Our interpreta- 
tion of this is that the absorption system at -190 km s -1 covers 
the quasar nuclear emission more completely than it covers 
the extended Ly a emission. We will return to this point in 
Section 3.3. 

The middle panel of Figure [3] shows that the Ly a emis- 
sion at velocity to +700 km s" 1 has only a marginally greater 
FWHM than the continuum fit level marked by a dashed line. 
This does not preclude a significant extended component at 
this velocity due to the fact that a minor extended component 
does not alter the FWHM greatly. For example in the imag- 
ing of G09, the z' band quasar FWHM is equal to the model 
PSF despite the fact that 13% of the flux is in the extended 
component. 

In order to determine the PSF we use the Gaussian fit re- 
sults described above, which are dominated by data outside 
of the Ly a emission line wavelength range. This makes the 
very large assumption that the only extended flux in the spec- 
trum is in the Ly a emission line and not in wavelengths cor- 
responding to quasar continuum. We will return to analyze 
this assumption in Section 3.3. Our method is very similar 
to the spectral PSF (SPSF) method of Moller (2000; see also 
Moller et al. 2000). The linear fits to the Gaussian centroids 
and FWHM show that these quantities are very slowly vary- 
ing functions of wavelength and therefore the values expected 
at each wavelength are well-determined from the linear fits. 

As a first test of PSF-subtraction we subtracted Gaussians 
with the width and centroid determined from the linear fit de- 
scribed above. The normalization of the Gaussians were set 
by the observed peak values. This test showed that there were 
substantial positive and negative residuals in the form of wig- 
gles with peaks up to ±4% of the Gaussian peak value. We 
checked standard star spectra from the same run and found 
these wiggles to be present in Gaussians subtracted off stel- 
lar spatial profiles and therefore these wiggles form part of 
the instrumental PSF. The wiggles had subtly different form 
for the standard stars compared to the quasar data (which was 
collected over several different nights) and therefore the stel- 
lar spectra could not be used to define the instrumental PSF. 
Instead we used the quasar spectrum at wavelengths equiv- 
alent to Lya velocities of +1150 to +2150kms _I . At this 
wavelength the quasar is not obviously extended, yet there 
is sufficient flux to determine a high S/N PSF correction to 
the Gaussian profile. Note that this procedure means that 
any true extended flux at this velocity range would be sub- 
tracted off across the full spectrum. We note that our two-step 
spectral PSF method is different to the "look-up table" ap- 
proach of Moller (2000). However, the results for the two 
methods will be comparable because both correctly account 
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FIG. 4. — Signal-to-noise ratio of extended flux as a function of wavelength 
(plotted as Ly a velocity) for the simple case of PSF-subtraction scaling the 
PSF to the peak value of each column. The S/N values have been smoothed 
by 5 pixel averaging. Regions with S/N < 10 in the extracted quasar spectrum 
are marked with grey shading as in Figures[2]and[3] Negative S/N residual on 
this plot indicates negative measured flux due to over-subtraction of the PSF. 
The most negative regions are those with low S/N in the quasar spectrum 
where the Gaussian fitting is not reliable. The only large velocity region with 
positive residual S/N is the « 1000 km s~' centred on the quasar systemic 
velocity. 

for the wavelength-dependent PSF centroid, the wavelength- 
dependent PSF width and shape and the pixelization of the 
data. Evidence for this will be seen in the cleanliness of the 
PSF-subtracted quasar continuum (similar to the very clean 
PSF-subtracted continuum in Weidinger et al. 2004 using the 
look-up table method). 

The final problem is one inherent to PSF-subtraction, how 
to scale the PSF in order to just remove the unresolved compo- 
nent without over-subtracting the extended emission. By sim- 
ply scaling the PSF peak to the observed peak, the center of 
the recovered extended emission inevitably has zero flux. The 
first step in this process is to determine at what wavelengths 
we have unequivocal evidence for extended emission. This 
was achieved by performing spectral PSF-subtraction with the 
PSF peak set equal to the observed peak such that there is zero 
flux in the extended component at the center. The residual flux 
from this process was then integrated along each column to 
determine the total extended flux per wavelength pixel. This 
total residual flux was then divided by the total noise in each 
column of the spectrum, taking into account sky and object 
contributions to the noise. This resulted in the S/N of ex- 
tended emission as a function of wavelength and is plotted in 
Figure [4] The only velocity region which shows positive ex- 
tended flux residual over a sizeable velocity range is between 
-400 and +470 km s" 1 . Even though the S/N of each column 
is only between 1 and 2, the integrated S/N across this veloc- 
ity range is 10 and highly significant. 

Recall that the velocity region +1150 to +2150kms _1 was 
used to determine the non-Gaussian component of the PSF, so 
any extended flux at that velocity would be subtracted off as 
part of the PSF. Whilst we cannot conclude that there is no ex- 
tended emission at > +500 km s" 1 (see also Section 3.3), we 
can say that there is a greater amount of extended emission 
located close to the quasar Ly a center than at longer wave- 
lengths. 

Therefore the PSF is scaled to the observed peak outside 
of the velocity region -400 and +470kms _1 . Within this re- 
gion the PSF scaling is determined using the observed peak, 
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FIG. 5. — Upper: One-dimensional extracted quasar spectrum (black) com- 
pared to the PSF-subtracted residual spectrum of the extended emission (red). 
Middle: Original 2D spectrum. Lower: PSF-subtracted 2D spectrum show- 
ing the extended Ly a halo. Positive spatial offset corresponds to the east and 
negative to the west. 

but with a correction if necessary to prevent the total flux of 
the extended component within the central 2 or 3 spatial pix- 
els being negative. Note that this is still a maximal subtrac- 
tion and it is possible that the extended component contains a 
greater fraction of the total flux. 

3.3. Results 

The results of the spectral PSF-subtraction process are 
shown in Figure[5] The upper panel compares the quasar spec- 
trum with the spectrum of the extended component. The mid- 
dle panel shows the original 2D spectral image and the lower 
panel the 2D spectral image extended residuals after PSF- 
subtraction. There are very few positive or negative residuals 
at velocities < -400 or > +500 km s -1 which shows how well 
our spectral PSF-modeling works. The extended emission ap- 
pears clumpy both spatially and kinematically. This structure 
will be investigated further below. 

One of the surprising things about the residual extended 
emission is that its flux is so low. Recall that in Section 2 
we showed that the extended emission identified in Subaru 
imaging contributes 13% of the total z' band flux, assuming 
a maximal PSF-subtraction. To compare to the extended flux 
in the spectrum we need to consider differential slit losses for 
the PSF and extended emission and the wavelength range over 
which to integrate. Differential slit losses for the ESI seeing 



*g 0.5 
5 0.4 
? 0.3 



* 0.5 



: Total 








|FW 


M:640 km s"' |: 









West 








jhJ 


|FWHM:770 km s"' 1 

, } 


East 










|FWHM:570 km s"' |: 







-2000 -1000 1000 2000 

Velocity (km s" 1 ) 



FIG. 6. — Ly a halo flux velocity profile for the eastern component (lower), 
the western component (middle) and the total (upper). Wavelength regions 
with S/N < 10 in the extracted ID quasar spectrum are marked with grey shad- 
ing as in Figures [2] and [3] The red curves show Gaussian fits to the velocity 
profiles. None of the velocity profiles are well fit by a Gaussian. 

based on the Subaru imaging show that we expect 10% of the 
total observed spectroscopic z' band flux to come from ex- 
tended emission. 

Direct measurement of the flux between 9000 and 9040 A in 
the upper panel of Figure [5] shows that the extended emission 
contributes 10% of the total flux over this narrow wavelength 
range containing the extended Ly a emission. Over a wider 
wavelength range of 9000 to 9 1 00 A, the extended component 
is only 5% of the total flux. Based on the Suprime-Cam z' 
band response curve, 39% of the flux from the quasar in this 
filter comes from 9000 to 9 100 A. Therefore, only 2% of the z' 
band flux in the ESI slit comes from the extended component, 
compared to the expected 10%. The majority of the extended 
emission observed in Subaru imaging is unaccounted for in 
our spectroscopy. 

Possible reasons for this missing extended flux are that the 
broad-line region/continuum flux (which we used to deter- 
mine the PSF) is also spatially resolved or that we have sub- 
stantially overestimated the PSF scaling for the narrow part 
of the Lya line. G09 proposed that 40% of the extended 
z' band emission was continuum flux from an extreme star- 
burst (M1450 = 23.9) in the host galaxy, based on their tenta- 
tive detection of extended emission at z r band. We have no 
evidence for extended continuum based on Figure H] Also 
the PSF-fit appears to do an equally good job at accounting 
for all the observed ESI flux at velocities corresponding to 
the broad Ly a emission (velocities -2000 to -500 km s" 1 and 
+500 to +2000 km s" 1 ) as it does at longer wavelengths where 
the quasar continuum dominates. Even the extreme starburst 
advocated by G09 would be a minor component in the broad 
Ly a wavelength region, so this argues against significant ex- 
tended continuum. It is more likely that we have overesti- 
mated the PSF scaling for the Ly a line, since we have no 
independent constraint on this normalization. In this case up 
to half of the narrow Ly a emission could be in the extended 
component. The fact that our seeing is substantially worse 
than that of the imaging observations is another complicating 
factor. 
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FIG. 7. — Radial profile of surface brightness for the Lya halo of 
CFHQS J2329-0301 based on ESI spectroscopy. Different curves are plot- 
ted for the east and west components. Also plotted is the surface brightness 
profile from narrow-band imaging (filter NB921) of the z = 6.6 Ly a blob 
Himiko of Ouchi et al. (2009). The 2 a limit arrow shows the 2 a surface 
brightness limit for the CFHQS J2329-0301 data. All curves are plotted out 
to a distance of one data point beyond the last 2 a detection. 

Even though our observations leave some puzzle about the 
bulk of the extended emission, it is interesting to consider the 
spatial structure and kinematics of the resolved Ly a emission. 
Figure [6] shows the velocity profile for the whole Ly a halo 
and for the two spatially distinct components on the west and 
east sides. The overall velocity structure for the two sides is 
quite similar. The complex structure means that neither com- 
ponent is well fit by the Gaussians plotted in Figure [6] The 
FWHM of these Gaussians do however give a good indica- 
tion of the velocity range of the bulk of the Ly a emitting gas. 
Both contain a pair of narrow peaks just redward of zero ve- 
locity (more prominent in the east component). 

Both the east and west show a broad dip at -190 km s" 1 due 
to partial covering by the z aDS = 6.412 Ly a absorption system 
discussed in Section 3.2. In the eastern component this dip 
reaches to zero flux at - 1 50 km s" 1 . The eastern component is 
more strongly absorbed than the west, consistent with the spa- 
tial offset marked in the lower panels of Figures [2] and [3] To 
constrain further the difference in spatial coverage of the ab- 
sorber, we estimated the velocity spread 8v of the absorption 
in the quasar and in the total extended emission spectra. At 
an optical depth r = 0.1, Sv = 500 km s -1 towards the quasar, 
whereas it is only 360 km s" 1 for the halo, the latter not ex- 
tending as far in the blue as that towards the nuclear emission. 
Using the four component fit of the Ly a absorption profile, 
we find a velocity spread 5v = 435 km s -1 for the two outer- 
most components. The total rest equivalent widths of Ly a 
in the quasar and extended spectra also roughly constrain the 
variations in the amount of gas in the z a b s = 6.412 absorber. 
We find w r = 1.31 ±0.04 and 0.72 ±0.07 A for the quasar 
and halo, respectively. This suggests a significant difference 
in the amount of absorbing gas towards the nucleus and the 
extended halo. 

Figure [7] shows the surface brightness profile along the slit 
in each direction from the quasar. Emission on both sides of 
the slit is visible up to 1.5 arcsec from the quasar, similar to 
the physical extent visible in Subaru imaging (Figure[T|i. The 
two components have similar surface brightness distributions, 
although the western component appears to have a sharp cut- 



off at distance 1.3 arcsec. 

Also plotted on Figure |7]is the radial profile of the z = 6.6 
Lya blob Himiko from the narrow -band imaging of Ouchi 
et al. (2009). The size and surface brightness profiles 
of these two objects are very similar, although the halo of 
CFHQS J2329-0301 is a factor of 2 brighter than Himiko. 
This will be discussed further in Section 4. 

4. DISCUSSION 

The Lya halo around CFHQS J2329-0301 has a physical 
extent of at least 15 kpc. By comparison, typical star-forming 
galaxies at z = 6 have a half-light radius of 1 kpc (Bouwens et 
al. 2004). The Ly a halo contains gas with a velocity range 
of » 1000 km s" 1 . The total emission line flux within the 1 
arcsec ESI slit is 8 x 10~ 20 W m~ 2 which gives a line luminos- 
ity of 4 x 10 36 W. Accounting for the fact that only half the 
flux identified in the Subaru imaging of Section 2 is located 
within the ESI slit and that the majority of the extended emis- 
sion from imaging is not identified in our spectroscopy, the 
total flux and luminosity are at least twice as great as this, i.e. 
>8xl0 36 W. 

The halo line luminosity is comparable to the most lumi- 
nous Ly a blobs known (Steidel et al. 2000). Christensen et 
al. (2006) showed that radio-quiet quasars at 3 < z < 4 have 
Ly a halos with luminosities typically 0.5% of the broad Ly a 
luminosity. By comparison, for radio-loud quasars (Heckman 
et al. 1991), this ratio is typically 5%. For CFHQS J2329- 
0301 the ratio is 5%, despite the fact it is not a strong radio 
source (Wang et al. 2008). 

Interpreting the physical nature of Lya halos is notori- 
ously difficult. Ly a is resonantly scattered so the observed 
spatial and kinematical distributions are not those where the 
photons were emitted. Ly a photons can be generated from 
several different mechanisms including shock heating of in- 
falling/outflowing gas or photo-ionization by stars or an AGN. 
For the case of CFHQS J2329-0301, we know that the center 
of the galaxy harbors a 2 x 10 M Q black hole accreting at the 
Eddington limit (Willott et al. 2010) which is generating a lot 
of Ly a emission within the quasar broad-line region (upper 
panel of Figure 2). The influence of the quasar UV emission 
extends 2.5 Mpc along our line-of-sight as evidenced by the 
transmitted flux at -2000 km s -1 in the upper panel of Figure 

2 (see also Willott et al. 2007). This makes the AGN explana- 
tion to power the Ly a halo a distinct possibility. The size and 
surface brightness of the observed halo is within the range 
predicted by Haiman & Rees (2001) for cold gas photoion- 
ized by a quasar. This implies the galaxy is at an early stage 
in its evolution due to the copious amount of cold gas located 
at large distances from the nucleus. The higher ratio of halo 
to broad-line luminosity compared to radio-quiet quasars at 

3 < z < 4 could be due to a greater covering factor of cold gas 
at higher redshifts, but a larger sample at z ~ 6 is required. 

What clues come from the gas kinematics? The gas is ap- 
proximately centred at the quasar systemic velocity when ac- 
counting for the absorption by the system at - 1 90 km s" 1 . The 
infall models of Villar-Martin et al. (2007) and Weidinger et 
al. (2004) both expect a spatial offset of the peak emission, 
a velocity offset up to 500 km s -1 and a velocity FWHM that 
varies systematically across the object. None of these features 
are observed for this Ly a halo. Therefore there is no evidence 
that the gas is dominated by either outflowing or infalling gas. 

The velocity FWHM of 640 km s" 1 reflects the velocity 
range of the gas where the Ly a photons are emitted and scat- 
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tered within the central 15kpc diameter of the host galaxy. 
This FWHM is much higher than that of the infalling Ly a 
gas in the z = 3 quasar Q 1205-30 (Weidinger et al. 2004) 
or of Himiko which has FWHM= 250 km s -1 . However, it is 
comparable to other Ly a blobs found at z = 3 which are sug- 
gested to be the sites of massive galaxy formation based on 
these large line widths (Matsuda et al. 2006) and to the quasar 
Ly a halos observed by Christensen et al. (2006). 

Making the possibly unreasonable assumption that the gas 
is in virial equilibrium at a typical distance of 7kpc with 
velocity dispersion a = 640/2.35 = 270kms _1 , one can de- 
termine a dynamical mass of M = 6 x 1O U M0 using M = 
(5/3)(3a 2 )R/G. Haiman & Rees (2001) show that even mas- 
sive halos such as that based on this dynamical mass can con- 
tain substantial quantities of cold gas in their inner regions 
due to the short cooling time. We note in passing that ellipti- 
cal galaxies in the local universe with a = 270 km s" 1 typically 
contain black holes with the same mass as CFHQSJ2329- 
0301 (Tremaine et al. 2002). This is in contrast to the results 
of Wang et al. (2010) which found SDSS z ~ 6 quasars have 
black hole mass to bulge mass ratios a factor of « 15 greater 
than locally. 

5. CONCLUSIONS 

We have spectroscopically identified a large and extremely 
luminous Lya halo around one of the most distant known 
quasars. The Lya properties do not show signatures of ei- 
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